Synthesis of small and large scale dynamos 



Kandaswamy Subramanian 

National Centre for Radio Astrophysics, TIFR, Poona University 
Campus, Ganeshkhind, Pune ^11 007, India. 

o. 

Abstract. Using a closure model for the evolution of magnetic corre- 
lations, we uncover an interesting plausible saturated state of the small- 
scale fluctuation dynamo (SSD) and a novel anology between quantum 
2^ | mechanical tunneling and the generation of large-scale fields. Large scale 

fields develop via the a-effect, but as magnetic helicity can only change 
CO | on a resistive timescale, the time it takes to organize the field into large 

scales increases with magnetic Reynolds number. This is very similar to 
the results which obtain from simulations using full MHD. 
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1. Fluctuating field dynamics 

I The dynamics of the fluctuating magnetic field B, is governed by the induc- 

tion equation. The velocity is assumed to be the sum of a Gaussian ran- 
dom, delta-correlated in time vy, and an ambipolar diffusion type component 
V£> = a[(V x B) x B)]. (Here a = r/(47rp), r is some response time, and p is 
the fluid density). Assuming that the magnetic field is also Gaussian random, 
Subramanian (1997, 1999; S99) derived closure equations for the longitudinal 
correlation function M(r, t) and the correlation function for magnetic helicity 
density, N(r,t). The random has a longitudinal correlation function T{r) 
and a correlation function for the kinetic helicity density, C(r). Defining the 
operators D(f) = (1/V 4 ) (d(r i f)/dr), and D(f) = (df /dr), we then have, 

M = 2D(t] T DM) + 2GM + 4aH; N = -2rj T H + aM, (1) 

where H = —DDN is the correlation function of the current helicity, G = 
—DDT is the effective induction. Also a = ao(r)+4aH(0, t) and ??t = ??+%( r ) + 
2aM(0, t) are functions resembling the usual a-effect and the total magnetic 
diffusivity. Here ao(r) = — 2[C(0) — C(r)] and rjo( r ) = T(Q) — T{r) and r] the 
microscopic diffusion. Note that at large scales r — > oo, a — > = — 2C(0) + 
4aH(0,t) and t]t — > f?oo = T) + T(Q) + 2aM(Q,t). The a-effect suppression is 
similar to the a-suppression formula first found by Pouquet et al. (1976); rjr is 
however enhanced by the growing field energy, as in ambipolar diffusion. 

1.1. Small-scale dynamo saturation 

The SSD problem, with C(r) = 0, a = has solutions with H{r, t) = 0, and was 
first solved by Kazantsev (1968). A transformation of the form \t(r) exp(2ri) = 
r 2 ^/rfrM(r,t), maps the problem of getting T > modes, into a bound state 
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problem in time-independent quantum mechanics (QM). For turbulent motions 
on a scale L, with a velocity v, bound states obtain provided the magnetic 
reynolds number (MRN) R m = vL/rj > R c m 60, and imply M growing at 
rate ~ v/L. Further, the bound-state eigen-function describes a field, which is 
strongly concentrated within the diffusive scale, r = ~ L{R m )~ 1 / 2 , and curved 
on the scale L. For a / 0, rj is simply replaced by an effective, time dependent 
diffusion rjo = rj + 2aM(0, t). So as the field (and M) grows, the effective MRN 
Ro{t) = vL/i]£>(t) is driven to the critical value R c . The final saturated state is 
obtained when R£>(t s ) = vL/(r) + 2<iMl(0, t s )) = R c ~ 60. So at saturation, the 
average energy density E B {t s ) = (3M L (0, i a )/8?r) = {3/2)(fn? /2)(L/vt)(1/Rc). 
For r ~ L/v, Eb is a small fraction ~ Rq 1 <C 1, of the equipartition value. If 
we intepret the saturated field configuration in terms of flux ropes with peak 
field Bp, thickness r^, and curved on scale L, Eb ~ (i?p/87r)Lr^/L 3 . Using 
r 2 d jL 2 -R^ 1 , and r ~ L/v, we then have B 2 /8ir ~ pv 2 /2, where, remarkably, 

the Rq 1 dependence has disappeared. So the SSD could saturate with the 
small-scale field of equipartition strength, being concentrated into flux ropes 

of thickness LR C 1//2 , and curved on scale L. 

1.2. Large scale dynamo as a tunneling problem 

For helical turbulence with C(0) 7^ 0, new generation terms arise at r 3> L, 
due to the a- effect, in the form M = .... + 4a 00 .ff and N = ... + OooM. These 
lead to the growth of large-scale correlations on a scale D, with a growth rate 
~ ctoo/D — r/oo/D 2 , as in the large-scale a 2 - dynamo. A special wavenumber, 
k p (t) = ojqo (t)/rjoo (t) , is also picked out for any quasi-stationary state. For 
such states, which obtain when one neglects slow resistive evolution (see below), 
(dN/dt) w 0, implying H w (a/2r] T )M. And if we define * = r 2 ^/rj^M, the 
equation (dM/dt) ps 0, can again be mapped to a QM potential problem, for 
the zero-energy eigen-state. However the modified potential now tends to a 
negative definite constant value of —a^/r]^ at large r and so allows tunneling 
(of the bound state) (see S99). In fact for r 3> L, one has an analytic solution 
M(r) oc r -3 / 2 J 3 / 2 (fepr), exactly as one would get if the large scale field, Bo, was 
random and force- free with V x Bo = A; p Bo. 

2. Helicity constraint and resistively limited growth 

The closure equations have also been solved numerically by Brandenburg and 
Subramanian (2000) (BS2K), adopting forms for T(r) and C(r) to match closely 
with the direct simulation of Brandenburg (2000) (B2000) (Run 5). One sees 
an initial exponential growth of the magnetic field, which terminates when its 
energy becomes comparable to the kinetic energy. Note that our closure equa- 
tions satisfy the helicity constraint N(0) = —2r/H(0). The numerical solutions 
show that, after some time t s , the current helicity (J • B) oc H(0,t), is driven 
to a constant value which however is such that |ci!oo| remains finite. A constant 
H (0, t) implies that the magnetic helicity (A • B) oc iV(0, t) grows linearly at 
a rate proportional to rj. During this phase the magnetic field correlations can 
extend to larger and larger scales. The corresponding magnetic energy spectra, 
Ey[(k, t) = (l/7r) Jq 00 M(r, t) (kr) 3 j\{kr) dk are shown in Figure 1. 
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Figure 1. Evolution of magnetic energy spectra. Note the propaga- 
tion of magnetic helicity and energy to progressively larger scales. The 
k~ 2 slope is given for orientation, (see BS2K) 



The resulting magnetic field is strongly helical (cf. section 1.2) and the mag- 
netic helicity spectra (not shown) satisfy |Am| ^ {2/k)Ey[. The development 
of a helicity wave travelling towards smaller and smaller k, as seen in Figure 1, 
is in agreement with the closure model of Pouquet et al. (1976). We have also 
checked that to a very good approximation the wavenumber of the peak is given 
by fcpeak(i) ~ k p (t), as expected from section 1.2, and it decreases with time be- 
cause c^oo tends to a finite limit and 7700 increases. Further, since the large scale 
field is helical, and since most of the magnetic energy is by now (after t = t s ) 
in the large scales, the magnetic energy is proportional to (B 2 ) ~ fc p (A • B), 
and can therefore only continue to grow at a resistively limited rate. These 
results are analogous to the full MHD case (B2000); the helicity constraint is 
independent of the nature of the feedback! In conclusion, our closure model 
with ambipolar diffusion type non-linearity provides a useful model, enabling 
progress to be made in understanding nonlinear dynamos. One now needs to 
think of ways, to relax the helicity constraint, (cf. Blackman and Field (2000), 
Kleeorin et al. (2000)), so that large-scale magnetic fields can grow fast enough. 
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